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new physics signals. Here we focus on data-driven backgrounds to a missing-energy-plus-jets 
search performed by CMS. Finally, we review recent progress in developing theoretical tools 
for high-multiplicity loop-computation within the BLACKHAT-library. In particular, we discuss 
methods for handling the color degrees of freedom in multi-jet predictions at NLO. 
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1. Introduction 



The production of a massive vector boson in association with jets is an important process 
at hadron colliders. The cross sections are large, the events are relatively clean, and they form 
significant backgrounds to many interesting physics signals. Because these processes are so well 
understood, both experimentally and theoretically, they are widely used to vaUdate or test new tools 
and methods. 

The last few years have seen continued progress in the perturbative description of high-multi- 
plicity processes. The next-to-leading order (NLO) QCD corrections for W + 4-jet production at 
hadron colliders were completed in 2011 [|], followed by the calculation of the same process 
with a W boson replaced by a Z boson At this conference we have shown preliminary re- 
sults for the W + 5-jet process. All these fixed-order QCD predictions have been obtained using 
of BlackHat [§ and SHERPA [|, g % There has also been a lot of progress in the com- 
putation of processes with a W and Z boson accompanied by jets at NLO accuracy, matched to a 
parton shower (see e.g. references in ref. [Q]). For such computations the high-multiplicity virtual 
matrix elements in QCD are a key ingredient, some of which have become available only recently. 
The production ofW + 2 jets has been computed using aMC@NLO [^. Virtual matrix elements 
provided by BlackHat have been used by different groups for such computations. The Sherpa 
implementation of the MC@NLO [ [Io| ] approach computes W boson production in association to up 



to three jets ||11|] using virtual matrix elements from refs. [|12[ [IJ. The POWHEG BOX [|14[ [15|] 



has been used recently to compute the Z + 2-jet process at NLO, matched with a parton shower 



[]16[], with the virtual matrix elements [|17[ [ISp also employed in ref. []19|]. 

NLO predictions improve the leading order (LO) results in various ways. NLO results show a 
reduced dependence on the unphysical renormalization and factorization scales, as compared to LO 
results. This improvement becomes more significant as the number of jets and, thus, the order in the 
strong coupling a,, increases. Further benefits include a better description of initial and final state 
radiation. The shapes of many kinematical distributions are better described at NLO. The precision 
offered by NLO calculations is needed for both the signal and background processes. 'Data-driven' 
methods very often rely on theory input for cross section ratios, in order to extrapolate backgrounds 
from a control region into a signal region for the same process, or to extrapolate from one process 
to a related process. NLO computations can also improve the theoretical precision for such ratios. 

We address here several new developments in fixed-order NLO computations of vector-boson 
production. First, we discuss the production of a massive vector boson in association with four 
jets at NLO as a signal at the LHC. Next, we discuss an application of NLO computations of 
vector-boson production to background estimation in searches for supersymmetry. Finally, we 
report on some technical developments that have been important ingredients in our latest multi-jet 
computations. 

* Speaker. 
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2. W + 4-jet and Z + 4-jet predictions 

We present results here for W + 4-jet and Z + 4-jet production, for which we have used 
BlackHat [Q] for both the virtual matrix elements and real-emission matrix elements. SHERPA 
^ 0, H] was used for the remaining pieces (Born and subtraction), as well as for the integration over 
the phase space. For both processes, the calculation of the real-emission tree amplitudes is ex- 



tremely challenging. In BlackHat, they are computed using on-shell recursion relations [|20p as 
well as compact analytic formulae given in refs. 



In figure |I] we display the transverse-momentum {pj) distribution of the first, second, third and 
fourth jets in Z + 4-jet events, with the decay to a lepton pair included. We have used Hj/2 as a cen- 
tral choice for both the factorization scale jXf and renormalization scale jiR, where Hj = Ej + Y^iPj 

and = ^Ml+{p'j''f. The sum runs over all paitons. A detailed list of the cuts, jet algorithm 
and parameters can be found in ref. To suppress the virtual-photon component, we apply a 
cut on the lepton invariant mass, requiring 66GeV < M^+g- < 116GeV. Although the interference 
between Z boson and photon exchange is very small with this cut, we include it for complete- 
ness. For the virtual part we use a leading-color approximation, which has been demonstrated in 
refs. [13, 19, ^ to be good to about three percent, as discussed further in section 4^ 



In the top pane of figure [T], the blue and black curves are the central LO and NLO predictions, 
respectively. In the middle pane the ratio is taken with respect to the central NLO result. Scale 
variation bands are displayed in orange for the LO and in gray for the NLO result. They are 
obtained as the envelope from varying }X = iIr = }Xfhy factors of 1/2, 1/ \/2, 1, y/l and 2 around 
the central scale value. The bottom pane displays the ratio of the Z + 4-jet distributions with respect 
to the + 4-jet ones, both at leading order and next-to-leading order. In contrast to the individual 
distributions, the ratios do not suffer from large NLO corrections. The shape of the Z/W^ ratio is 
explained by the dominance of the valence distributions, particularly u{x), at large x. The Z has a 
significant coupling to the u quark, while the couples only to the d quark, among the valence 
quarks. As the jet pj increases, parton-distribution functions (PDF's) of higher ;c are probed, where 
u{x)/d{x) rises. The Z/W^ ratio is flat because both the W and Z couplings are dominated by the 
coupling to the same u{x) PDF. These sorts of ratios are very useful in data-driven methods. 

The shape differs significantly between the central LO and NLO predictions for the pj distri- 
butions of the first, second and third jets. The scale variation is, as expected, much smaller at NLO 
than it is at LO. 



3. Z/ 7 ratios 

The production of jets in association with a Z boson that decays into a neutrino pair is an 
irreducible background for missing energy plus jets (METJ) searches. Missing energy signals arise 
from the production of new particles that escape the detector unobserved. A typical example is the 
lightest supersymmetric particle (LSP) in /^-parity conserving models of supersymmetry. As it is 
an in^educible background, the impact of the Z(— )• vv) + jets processes (METZJ) must be estimated 
cai^efully. Typically data-driven methods are used with only ratios provided by theory. There ai^e 



several strategies to estimate this background (see e.g. refs. [23, 241 23, E6]), using related processes 
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Figure 1: Distribution in pj for the first, second, third and fourth jets in Z + 4-jet events. 



in which the production of the Z boson and its subsequent decay into neutrinos is replaced by a 
more accessible signature: 

• Z (— )• //)+jets. The advantage of this process is that the production dynamics are very similar 
to the one in the neutrino process. However, the statistics is smaller than Z — )• vv by a 
factor of six for each lepton flavor, and is further reduced by the experimental acceptance 
and identification efficiency for each charged lepton. 

• W (— )• /v)+jets. This process has a cross section higher by a factor of six, but suffers from 
contamination from tt, and potentially, from the new physics one aims to measure. 

• 7+jets. Replacing the Z decay to neutrinos by a photon yields a cross section higher by a fac- 
tor of four to five, but the production dynamics are different and a reliable theory prediction 
is needed to obtain a solid conversion factor. 



Here we present NLO calculations of the ratios needed for the photon-based estimation strat- 



egy. The CMS collaboration has studied and used [^, W-boson and photon production in 
association with jets to estimate METZJ backgrounds [U, ^]. Photon production in association 
with jets has also been studied in ref. [25] and used by the ATLAS collaboration [ ^ ] in their 
data-driven estimates of the METZJ background. 

Partonic calculations involving photons develop infrared singularities when the photon is emit- 
ted coUinear to a quark. These divergences can be avoided by imposing a standard cone isolation 
for the photon. However, this strategy requires the use of photon fragmentation functions which 
are extracted from data and are not known very precisely. Another approach, proposed by Frixione 
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Figure 2: Ratio of the pj distribution of the first jet in Z+ 2-jet and 7+ 2-jet production. The dotted blue 
line represents the LO ratio, the solid black line the NLO result, while the ME+PS result is represented by 
the dashed red line. 



[28], removes the need for the fragmentation functions but is difficult to implement in an experi- 
mental measurement. For the estimation of the ratio of 7 + 2 jets to Z + 2 jets we use the Frixione 
isolation. We estimate the difference between this isolation and the CMS cone isolation using a 
code by Gordon and Vogelsang [ ^ ] and JetPhox ||3^, 31]. We find that the difference between the 
standard cone isolation and Frixione's isolation is small, and decreases as the transverse energy 
of the photon increases, which is the case of interest for this study. For more details we refer to 
ref. [||]. 

The estimation of uncertainties on our result for the ratios (as for most NLO predictions for 
ratios) is difficult. The commonly used method of varying the renormalization and factorization 
scales yields a very small estimate, because most of the scale-dependence is correlated between the 
numerator and denominator. A similar effect occurs for the PDF-eiTor propagation. To estimate 
the error, we computed the ratio using the matrix-element-matched-to-parton-shower (ME-i-PS) 
method of SHERPA. We took the difference between the NLO and ME-hPS ratios as an estimate 
of the uncertainty of our result. An example of such a ratio, for cuts relevant to a CMS search using 
the 2010 LHC data 124], is shown in Figure ^ as a function of the pj of the first jet. The NLO 
and ME-i-PS predictions for the Z to / ratio track each other well across the whole range of jet pj. 
Other observables exhibit similar behavior [|3^], leading to the conclusion that the photon-plus-jets 
process gives a good handle on the Z — ^ vv background. 



We have recently extended this computation [ [33| ] to the ratio of Z+ 3-jet to 7+ 3-jet pro- 
duction. This work was also instigated by the CMS collaboration [|^] to help them estimate the 
uncertainty in the Z — )• vv background to new physics searches, in particular for the more stringent 
cuts used in analyzing the 2011 LHC data [26]. We have computed the ratio of the two processes 
for different sets of cuts on two kinematic variables Hj^^ and MET^^^^. Hj^^ is defined as the sum 
of the transverse energy of jets with pj > 50 GeV and < 2.5, while MEP^^^ is defined as the 
modulus of the vectorial sum of the transverse momenta of the jets with pr > 30 GeV and |t] | < 5. 

The scale Q of the cuts on Hj^ and MET^^^^ in the analysis is large enough that one may woiTy 
about large logarithms of the form log(2/p^'"). One way of assessing whether these logarithms 
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Figure 3: Ratio of the Z + 3 to Z + 2 jets ratio as a function of //f " and //f " - MEV""^" . 



are large is to look at the ratio of Z + 3 jets over Z + 2 jets. A value too close to one could result 
from large logarithms spoiling the perturbative expansion. Figure ^ shows this ratio as a function 
of //j?'' and //j?" -M£■^j'='^ In regions where this ratio is close to one, we cannot be sure about 
the validity of the perturbative expansion. Such large effects may nonetheless largely cancel out 
in Z+jets/y+jets ratios, leaving those predictions reliable in spite of them. It is worth noting that 
the validity of the NLO prediction in these regions, where the perturbative expansion might be 
unreUable, can be tested experimentally by measuring such ratios in the 7+ jets samples. 

We find that the theoretical uncertainty on the conversion between photons and Z bosons is 
less than 10% for events with either two or three associated jets, allowing the photon channel to 
provide an excellent determination of the Standard Model missing-ET- + jets background. 



4. Color Automation 



Scattering processes at hadron colliders are dominated by the strong interactions. The previously- 
discussed processes of vector-boson production in association with jets are typical examples. For 
events with kinematics in the perturbative regime, cross sections and distributions can be predicted 
from first principles, allowing detailed comparisons between theory and experiment. To this end we 
focus on the precise perturbative description of the interactions of colored partons within quantum- 
field theory. 

The Lagrangian interaction terms of colored fields make perturbative computations very chal- 
lenging. A successful computational strategy disentangles the various dynamical degrees of free- 
dom, such as kinematics, spin and color quantum numbers. Targeted algorithms can then be de- 
vised to deal efficiently with the components of the computation. This approach is realized in 
many matrix-element generators and has been particularly useful for recent multi-jet computations 
at NLO. Alternatively, one can treat color simultaneously with the kinematic variables [Q]. This 
method, implemented as a Berends-Giele recursion [^], is similar to methods used in other auto- 
mated LO computations [36, 38, 3^. 

In BlackHat we separate color and kinematic degrees of freedom early on. We then compute 
color-ordered tree- and loop-level objects, so-called primitive amplitudes. These amplitudes are 
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assembled into the full matrix elements only at the end of the computation. The required color 
weights and interference matrix are precomputed. We use modern on-shell and unitarity techniques 
in order to compute the color-ordered building blocks: tree amplitudes and, subsequently, loop 
amplitudes. 

Focusing on these particular color-ordered objects has several benefits. First of all, these 
components have a simpler analytic structure than the full ampUtude; for example, they depend 
on a reduced set of kinematic invariants. This property implies that a smaller number of unitarity 
cuts needs to be computed in an on-shell approach. The simpler analytic structure also leads to 
improved numerical stability of the amplitudes. 

The second main benefit of the color-ordered approach is that it can be exploited to yield sig- 
nificant efficiency gains in the numerical phase-space integration [pj|]. If we consider the number 
of quark flavors «j to be of the same order as the number of colors Nc, and take the limit as both 
become large, then the color-summed (virtual) cross section can be expanded in powers of l/N^. 
The 1/A'^^^-suppressed terms are numerically quite small. Thus, for a fixed integration error, fewer 
evaluations are needed for these parts of the cross section. That is, for most phase-space points only 
a small subset of all color-ordered amplitudes must be computed, namely those that contribute to 
the leading-color term in the cross section. 

In this talk we discuss a method to automate the color-ordered approach in loop computations. 
The key question is: how we can express generic matrix elements as linear- combinations of color- 
ordered objects. We also discuss the quantitative impact of subleading-color terms on differential 
cross sections. In particular, we consider the distribution of the fourth jet pj in the state-of-the-art 
NLO results for W + 4-jet production [|13|, [T], g]. 



4.1 Partial amplitudes from primitive amplitudes. 

The standard "trace-based" color decomposition of a one-loop QCD amplitude is into a set 
of color structures involving traces or open strings of Nc x Nc SU (Nc) generator matrices T"*, one 
for each external gluon. The open strings terminate on fundamental indices coiTcsponding to the 
external quarks and antiquarks in the process. The coefficients of these color structures are called 
partial amplitudes. The partial amplitudes are in turn built from color-ordered primitive amplitudes, 
but the precise relations can be laborious to determine in the general case. The aim of this section 
is to describe an algorithm [Q] for determining these relations. 

For many specific processes, the decomposition of the one-loop matrix elements and the color- 
summed virtual cross sections in terms of primitive amplitudes is well known. The explicit color 



decompositions of all-n matrix elements of purely gluonic processes can be found in refs. [ 140| , |41[ ] 
and the ones with a single of quaik line in ref. [p2|]. Furthermore, the decomposition of the four- 
quark amplitude was given in refs. [|^, and that of the four-quark process with an additional 



gluon in [ ]44| , [45[ ]. The decomposition of QCD scattering amplitudes with six and seven partons, 
including either four or six quarks, have been given only recently, using the algorithm reviewed 
here [g]. 

The algorithm is based on analyzing Feynman diagrams and their inherent color information. 
We consider a specific subprocess in QCD. We generate the full set of Feynman diagrams of the 
colored loop amplitude. Next, we associate a linear combination of Feynman diagrams to partial 
amplitudes. To this end, we dress Feynman diagrams with color matrices and sum over repeated 
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(internal) color indices. Partial amplitudes are defined as coefficients of particular products of color 
matrices. Feynman diagrams ai^e associated with a given partial amplitude if they contribute to its 
defining color structure. Simple factors (powers of Nc, signs and integers) are generated by the 
color-index summation, and these factors enter the relation between each Feynman diagram and 
each partial amplitude to which it contributes. 

In a second step, we associate a linear combination of Feynman diagrams with primitive am- 
pUtudes. To do this, we again dress Feynman diagrams with color matrices; however, this time we 
associate adjoint representation color charges to quarks as well as gluons. Again we sum over the 
internal color indices. We define the primitive amplitudes as the leading-color single-trace partial 



amplitudes [40] for a given cyclic ordering of the external gluons and (adjoint) quarks. A primitive 
amplitude is thus associated with the Feynman diagrams that contribute to a particular single-trace 
color structure. For further details, and refinements when deaUng with quark lines, we refer the 



reader to the original literature [40, 41, \2 



Once we have the explicit expressions of both the partial amplitudes and the primitive ampli- 
tudes in terms of linear combinations of Feynman diagrams, we can manipulate these two sets of 
equations. We express the partial amplitudes in terms of the primitive amplitudes by solving the 
linear set of equations to eliminate the explicit dependence on the individual Feynman diagrams. 
In this step, we observe a redundancy of the linear equations we have to solve. This redundancy 
implies that certain linear- combinations of primitive amplitudes add up to zero. That is, we find 
non-trivial relations between primitive amplitudes. Thus, we may find a set of equivalent color 
decompositions of a given partial amplitude. 

We point out some key features of the implementation of the above algorithm. Because it is 
based on Feynman diagrams, the algorithm is limited in parton multiplicity due to the rapid growth 
in the number of diagrams. (However, the algorithm only has to be performed once and for all for a 
given process, not each time an amplitude is evaluated.) We find that there is no serious obstruction 
to carrying out the decomposition for processes with up to eight external partons. Processes with 
zero or two external quarks are computationally the most challenging; however, explicit formulas 
are known for these cases. In addition, it turns out that many diagrams (e.g. diagrams involving 
four-gluon interactions) may be dropped from the start, significantly reducing the computational 
load. For an efficient implementation, we find it convenient to use QGRAF [HH] for diagram 



generation and the computer algebra package FORM [47] for summing over internal color indices. 

We conclude this section with a brief discussion of the relations between primitive ampli- 
tudes, which appear as a byproduct of the above algorithm. At low multiplicity, one can identify 
the anti-symmetry of the color-ordered gluon-quark-quark vertex as the origin of the relations be- 
tween primitive amplitudes. Explicit examples of this may be found in ref. [§]. The relations 
obtained between primitive amplitudes can be used to optimize caching, and thus, how long the 
subleading-color contributions take to evaluate numerically. Finally, it seems likely that under- 
standing the relations between primitive amplitudes will help to establish all-n formulae for the 
color decompositions. We have not explored this direction further. 

4.2 Quantitative impact of subleading-color terms. 

We now discuss the quantitative impact of subleading-color terms on NLO predictions. We 
focus on the distribution in the fourth jet pj- in W + 4-jet production at the LHC, since it is a key 
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Figure 4: A comparison of the full and leading-color virtual contributions to the pj distribution of the fourth 
jet in W"+4-jet production at the 7 TeV LHC. 



observable. Predictions of jet pj- -distributions have already been given in ref. [jl|]. In that work a 
leading-color approximation was used for the virtual parts; the remaining real, Born and subtraction 
terms were computed to all orders in the color expansion. Here we show results including the 
remaining subleading-color corrections The systematics we observe match our earlier results 
on subleading-color contributions to W + 3-jet production [ 13]. 

We use the same basic setup of matrix-element generators, Monte Carlo integration, scale 
choices and cuts as discussed in section |2[ There is no unique definition of the leading-color terms; 
definitions may differ by reassigning subleading-color terms. A detailed discussion of the leading- 
color approximation we use here is provided in ref. [^]. 

In figure ^ we compare the full-color versus leading-color virtual contributions to the pj dis- 
tribution of the fourth jet in W^-i-4-jet production. Also displayed are the subleading-color contri- 
butions by themselves, labeled as 'full-minus-leading-color'. In order to obtain the full parton level 
differential cross section one must add the real and Born contributions in the usual way. The size 
of these contributions have been given recently [l3|] ; the leading-color virtual part accounts for 
about 20% of the cross section. 

Figure ^ confirms that the subleading-color contribution is suppressed uniformly by almost a 
factor of 10 over the displayed range of pj- The suppression appears consistent with the expected 
factor of y/N'^ with Nc = 3. With the leading-color virtual part accounting for about 20% of the 
total cross section, the subleading-color virtual part amounts to less than a 3% correction to the 
total. 

A possible exception to the uniform suppression would be near a zero of the leading-color 
virtual cross section. Such zeros are not excluded on general grounds, but they would have to 
survive the sum over a large number of different helicity configurations. A priori, there is no 
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reason to assume that the vanishing of the leading-color contribution forces also the vanishing of the 
subleading-color contribution. Close to a putative zero of the leading-color contribution we would 
expect a relative enhancement of the subleading-color piece in the virtual matrix elements. Even if 
a zero were to appear in a special kinematic configuration, the subleading-color virtual terms are 
expected to keep their relative size with respect to the more inclusive total cross section and typical 
differential cross sections. Certainly figure ^ does not show any evidence for suppression of the 
leading-color cross section. 

Although the leading-color results could point to kinematic configurations where subleading- 
color contributions might dominate the virtual matrix elements, only their explicit knowledge al- 
lows to determine their impact on the full cross section. Similarly, the size and uniformity of the 
subleading-color terms justifies the use of leading-color approximations for many multi-jet observ- 
ables. Nevertheless, explicit computations are of great importance for a definitive understanding of 
such multi-jet observables. 
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